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Abstract 

The jet structure in events with Drell-Yan-produced W bosons is discussed, and a 
new model for describing such event is presented. The model is shown to explain recent 
measurements of the T4^-jet rapidity correlation and predicts a transverse energy flow 
at high W rapidities (corresponding to probing small-x partons in one of the incoming 
protons) higher than conventional parton cascade event generators. 
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1 Introduction 



The description of small-x partons within hadrons has attracted a great deal of interest, 
especially after the measurements of the proton structure function F2 at x values down to 
10~^ at HERA, where a substantial increase was found ||l|,^. Although the small-x rise was 
first predicted by the so-called BFKL evolution equation [^,^ it soon turned out [^,^ that 
it could also be explained in terms of the conventional Altarelli-Parisi (DGLAP) evolution 
equations [^HlC|] . Instead, much of the focus has been directed towards the study of hadronic 
final states in deep inelastic lepton-hadron scattering (DIS) events at small x, and it has 
been suggested that the large flow of transverse energy in the proton direction found in 



such events is a signal of BFKL dynamics |11| 



Much can also be learned from comparing data with different models implemented in Monte 
Carlo event generators. So far it has been shown that generators built around a conventional 
DGLAP-inspired initial-state parton showers, such as Pythia 13| and Lepto |14|, with 



strong ordering in virtuality, completely fail to describe things like the forward transverse 



energy flow at small x, while a generator such as Ariadne |15| - although not implementing 
BFKL evolution, but sharing with it the feature that emissions are unordered in transverse 
momenta - describes such event features quite well 1^. 



Besides deep inelastic lepton-hadron scattering, Drell-Yan production in hadron-hadron 
collisions is one of the cleanest probes of hadronic structure. Recent results from the 



DO collaboration |18| at the Tevatron shows a surprising feature of events with Drell-Yan- 
produced W bosons, namely the decorrelation in rapidity between the W and the associated 
jets. Although the typical x-values probed in W events at the Tevatron is on the order of 
yJm\iflS 80/1800 0.04, for large rapidities of the W , one of the incoming partons has 
a much smaller momentum fraction of the proton (e.g. yw ~ 2 gives xi ~ 6 x 10~^ and 
xi ~ 0.3). Therefore it could be worth while to take the experience gained from studying 
small- 2; final states at HERA and try to apply it to large rapidity PF-production at the 
Tevatron. 

In this paper the Dipole Cascade Model (DCM) [^,|2^, on which the Ariadne program is 
built, is extended to also model the jet structure of Drell-Yan production events. The main 



feature of the DCM for DIS [21| is that gluon emission is treated as final state radiation 
from the colour dipole formed between the struck quark and the proton remnant as in 
fig. |l]a. In this way there is no explicit initial state radiation, and the proton structure 
enters only in the way the dipole radiation is suppressed due to the spatial extension of the 
proton remnant. This approach has some problems when it comes to describing features 
particular to the initial state, such as the initial splitting of a gluon into a q^^-pair. 

The simplest extension of the DCM to also treat Drell-Yan production would be to treat 
gluon emission as final-state radiation from the colour dipole formed between the two 
remnants, as in fig. [l|b. However, as is seen from fig. |2|, the leading order W^+jet diagrams all 
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Figure 2: The leading- order Feynman diagrams contributing to W + jet production corre- 
sponding to the annihilation (a) and Compton (h) diagrams. 

correspond to initial-state radiation (except for the last one, which is the least important). 
And in particular it is clear that if the gluon emission is treated as final-state radiation 
between the two remnants, it would be difficult to explain the contribution to the transverse 
momentum of the W from the diagrams in fig. |2[ 

In previous work the DCM was extended to also include the boson-gluon fusion diagram 
in DIS, which can be viewed as a special case of initial-state gluon splitting into qq. In 
this paper, this approach is further developed into a general inclusion of initial-state gluon 
splitting into qq. Also, a way is presented of taking into account the contribution to the 
transverse momentum of the W from the gluon emission, formulated it terms of radiation 
from the colour dipole between the two proton remnants. 

The layout of this paper is as follows. In sections § and |3| the treatment of gluon emission 
and initial state gluon splitting is presented. In section results for the M^-jet rapidity 
correlation from the improved DCM is compared to a leading-order calculation and with the 
conventional DGLAP-inspired initial-state parton shower approach of Pythia. Also some 
predictions are given for the transverse energy flow in high rapidity W events at Tevatron 
energies. Finally, in section the conclusions are presented. 
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2 Gluon Emission 



The DCM for e'^e'" annihilation and deep inelastic lepton-hadron scattering is described 
in detail in refs. [19-21| and only a brief summary of the features important for this paper 
will be given here. 



The emission of a gluon gi from a qq pair created in an e~^e~ annihilation event can be 
described as radiation from the colour dipole between the q and q. A subsequent emission 
of a softer gluon g2 can be described as radiation from two independent colour dipoles, one 
between the q and gi and one between gi and q. Further gluon emissions are given by three 
independent dipoles etc. 

In DIS, the gluon emission comes from the dipole stretched between the quark, struck by the 
electro-weak probe, and the hadron remnant. The situation is the same as in e^e~ above, 
except that, while q and q are both point-like in the case of e^e~ , the hadron remnant in 
DIS is an extended object. In an antenna of size I, radiation with wavelengths A <C I are 
strongly suppressed. In the DCM, this is taken into account by only letting a fraction 

a = fj./k± (1) 

of the hadron remnant take part in the emission of a gluon with transverse momentum k± , 
where is a parameter corresponding to the inverse (transverse) size of the hadron. 



The phase space available in dipole emission is conveniently pictured by the inside of a 
triangle in the K-y plane, where k = In A;^ and y is the rapidity of the emitted gluon as in 
fig. 1^. In these variables the dipole emission cross section also takes a particularly simple 
approximate form: 

da oc asdndy. (2) 

In DIS, assuming that the hadron is coming in with momentum (using light-cone coordi- 
nates) {P+, 0, 0), and is probed by a virtual photon (— 0), the triangular area comes 
from the trivial requirement 

k+g = k±ey < P+ 

k_g = k±e-y < (3) 

The condition that only a fraction of the remnant participates in an emission means that 

k+g < (m/A:x)P+ (4) 

and translates into an extra cutoff in the phase space corresponding to the thick line in 
fig. ^. This should be compared to the initial-state parton shower scenario, where gluon 
emission is given by 

_ 2asl + fg{x/z) dz dQ^ 
"^""^ " Svr l-z fg{x) z Q2 • 
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Figure 3: The phase space available for gluon emission in DIS (thin lines) and the extra 
restriction due to the extendedness of the proton remnant (thick line). The dotted line 
corresponds to a line of equal suppression due to the ratio of parton density functions 
entering into a conventional parton shower scenario. 

Identifying the ratio of structure functions in eq. (^) (corresponding to the dotted hne 
of equal suppression in fig. ^) with the extra cutoff (^ in the DCM, the two models are 
equivalent in the low-fc^ limit. 

As mentioned in the introduction, the simplest way of extending the DCM to describe 
gluon emissions in Drell-Yan events is to describe it as radiation from the colour dipole 
between the two hadron remnants. One problem with this approach is what to do with 
the transverse recoil from the gluon emission. In e~^e~, this recoil is shared by the q and 
q. In DIS, since only a fraction of the remnant is taking part in the emission, only that 
fraction is given a transverse recoil, resulting in an extra, so-called recoil gluon The 
corresponding procedure for Drell-Yan would be to introduce two recoil gluons, one for 
each remnant. However in that way it is impossible to reproduce the transverse momentum 
of the W as given by the 0{as) matrix element. 



The 0{as) matrix element for q + q ^ W + g production takes the form |23 



tu 



where s, t and u are the ordinary Mandelstam variables satisfying s + t + u = m^,. In order 
to reproduce this in a parton shower scenario, where the gluon is emitted "after" the W is 
produced, this has to be convoluted with the parton density functions, and the lowest-order 
VF-production matrix element, again convoluted with the relevant parton densities, has to 
be factored out. This introduces some ambiguities, which are solved by assuming that the 
rapidity of the W is the same before and after the gluon emission, resulting in the following 
cross section, expressed in the transverse momentum k'j_ and rapidity yg of the emitted 
gluon 

dag _ 2as fqjx'g) fqjx'g) 



dygdk]_ 37r fg{Xq) fq{Xq 



X 
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Figure 4: The phase space available for gluon emission in W production (thin lines) and the 
extra restriction due to the extendedness of the proton remnants (thick lines). The shaded 
triangle corresponds to the phase space area covered by the W . 



{k]_ + k±m±e^y){kl_ + k±m±e~^y){kj_ + m5_ + k±m±{e^y + e'^y)) ' 

where Ay = Vg — Hw-, w-l = k\ + m'yy, yw the rapidity of the W and Xi and x[ the energy- 
momentum fractions carried by the incoming partons before and after the gluon emission 
so that 

_ mwey^ / _ m±ey^ + /c^e^f 

mwe-y^ , m±e~y'^ + k^e'y^ ^ ^ 

^q = 7^ ' (9) 



assuming the q coming in along the positive z-axis and a total invariant mass of V^- 



In the limit k"^ ^ "i^/) (0) reduces to the simple dipole emission cross section in eq. (^), 
so it is clear that the strategy outlined above is a good leading log approximation. It also 
turns out that it is fairly simple to correct the first gluon emission so that, disregarding the 
ratios of parton densities, eq. (0) is reproduced. 

The ratio of parton densities in eq. (0) is instead approximated by the suppression of the 
phase space introduced for DIS in eq. (^), which in this case corresponds to suppressions 
on both sides of the triangle, as in fig. ^. One problem with this procedure is that the k'j_ 
of the gluon and hence of the W is limited by this suppression to 



kl < fi^fs/l, (10) 

which, with /i ~ 1 GeV, gives k± <, 30 GeV. To be able to describe high-A;_|_ W production, 
it is clear that the sharp cutoff in fig. ^, which in any case is an oversimplification, must be 
replaced by a smooth suppression. In |21| it was shown that introducing a power suppression 
in the disallowed regions in fig. ^ does not influence the general event shapes in DIS; it is 
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Figure 5: The transverse momentum spectrum of the W at the Tevatron. The full line is the 
prediction of the 0{as) W + g matrix element as implemented in Pythia using CTEQ2L 
parton density functions. The matrix element calculation was cut off at k±w = 10 GeV to 
avoid divergences. The dashed and dotted lines are the predictions of the gluon emission in 
the DCM with (3 = oo and (3 = 2, respectively. 



clear, however, that such a power suppression would greatly influence the high-fc^ spectrum 
of the W. 



First, however, the way of obtaining a transverse momentum of the W in the gluon emis- 
sions must be formalized. It is clear that, in the first emission, the gluon corresponds 
unambiguously to initial state radiation, and hence the k± of the gluon must be balanced 
by the k±_ of the W . In further emissions this is not the the dipole radiation is 

a coherent sum of the emission from the incoming partons and the outgoing, previously 
radiated, gluon. It is therefore argued that only gluon radiation that takes place close to 
the W in phase space should be able to influence the k± of the W:, only gluons emitted 
in the shaded region of fig. |^, corresponding to Pg+ < Pw+ and Pg- < Pw-i will have 
their transverse recoil absorbed by the W. Outside this region the transverse recoil will be 
treated as in the DIS case above. 



Figure ^ shows the 0{as) W + g matrix element prediction (as implemented in Pythia) 
of the k±_ spectrum of the W at the Tevatron, compared to the modified DCM described 
above. With a sharp cutoff in the phase space, it is clear that the DCM cannot describe the 
high-fex tail of the spectrum. Instead, a smooth suppression can be introduced by changing 
eq. (|^), allowing a larger fraction a' of the remnant to take part in the emission with the 
probability 

corresponding to a smoothening of the theta function suppression in fig. ^, giving a power- 
suppressed tail. As seen in fig. |5|, using (3 = 2 describes well the high-fej^ spectrum obtained 
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from the leading-order calculation using the CTEQ2L[| ||2^ structure function parametriza- 
tion. In the following, this value of (3 will be used, unless stated otherwise. 

This concludes the description of the W + gluon jet in the DCM. However, at small x, the 
gluon density in the proton becomes very large, and the Compton diagrams in fig. |2|b are 
dominating. 



3 The Compton Diagrams 



The matrix element for the Compton diagrams looks like |2|] 

st 

The s-channel diagram is, of course heavily suppressed for small k±, and looking only at 
the t-channel diagram, convoluting with parton densities and factoring out the zeroth order 
W production cross section, as in the gluon emission case above, the normal leading log 
initial-state parton shower cross section for the splitting of an incoming gluon into a qq 



pair |25[] is obtained: 

4^ fqix) z Q 



da, = ^iz^ + il-zf)l4^^^, (13) 



where = —t and z = rn^/s. 

In the DCM, however, there is no initial state gluon splitting into qq, and, just as in the 
case of final-state g ^ qq splitting this process has to be added by hand to the DCM. 



The simplest way is to introduce the initial-state g ^ qq splitting in the same way as in |26], 
as a process competing with the DCM gluon emission described above. The competition is 
as usual governed by the Sudakov form factor using ordering in k\. Rewriting eq. ( p!3| ) in 
terms of the transverse momentum k\ and rapidity of the outgoing quark, the probability 
of the first emission to be an initial-state g ^ qq splitting at a certain k^ and Uq is given 
by 

dPq{k\,yq) ^ daq{kj_,yq) ^ 
dk\dyq dk\dyq 



/fcL^ax 2/ f ^ / ( daq{kl,y'q) dag{kl,y' , ^ 

where the second factor is the Sudakov form factor, corresponding to the probability not 
to have any emission of gluons or gluon splittings above the scale k\ . 

^The CTEQ2L structure function parametrization is used in all analyses in this paper where applicable. 
None of the conclusions in this paper were found to be sensitive to this choice. 
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Technically, the extra process is implemented as follows. If the quark going into the hard 
interaction on one side is a sea-quark, the remnant on that side is allowed to "radiate" the 
corresponding antiquark according to eq. (13). After such an emission, the remnant is split 
in two parts according to the prescription described in ref. [21|, one of which forms a dipole 
with the "emitted" antiquark while the other retains the dipole colour connection of the 
original remnant. 



If the first emission is a g ^ qq splitting, the full 0{as) matrix element is used, and the 
rapidity of the W is assumed to be the same before and after the emission, as in the gluon 
emission case. A splitting later on in the cascade, the kinematic is fixed by requiring the 
non-radiating remnant to be unchanged. In all cases, the transverse momentum of the 
struck system will of course balance the k± of the emitted antiquark. Note that only one 
initial state g ^ qq splitting is allowed per remnant. This is a good approximation, since 
a second such splitting is heavily suppressed by the parton density functions. 

This procedure can be used not only in the case of W production, but for all processes with 
a hadron remnant present. In particular it can be (and i^ applied in the DIS case. 

In fig. ^a, the W k± spectrum at the Tevatron is shown, using the full 0{as) matrix element 
(as implemented in Pythia) and using the modified DCM with the initial state g ^ qq 
splitting as implemented in Ariadne^. Clearly the DCM does a good job of reproducing 
the high-A;^ tail of the distribution. In fig. ^b the DCM is compared with the two parton 
shower approaches of Pythia, one using only parton showers and one using first-order 
matrix elements with parton shower added. Since the DCM is a leading-log cascade, except 
that the first emission is uses the full matrix element, it smoothly interpolates between the 
pure parton shower description, which should be a good approximation for small k±, and 
the matrix element description, which is good for high k±, but has to be cut off at small 
k± to avoid divergences. 



4 Results and Predictions 



In ref. it was found that, when looking at the rapidity of the balancing jet in high-A;_L 
W events at the Tevatron, no correlation with the W rapidity was found, while a leading 
order and a next to leading order calculation predicted a strong correlation. It was also 
found that a preliminary implementation of the DCM model described here reproduced 
data fairly well and only gave a very weak correlation^. 

^This is the default in Ariadne version 4.06 and later. 

^AU results labelled Ariadne or DCM are actually generated using the zeroth-order W production 
in Pythia, with the CDM added and, where indicated, using the string fragmentation implemented in 
Jetset p2| . 

^The results presented here differ from the ones in ref. due to a bug introduced in the initial-state 
g qq splitting in the preliminary version of Ariadne used in that paper. 
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k±w (GeV) 



Figure 6: The transverse momentum spectrum of the W at the Tevatron. In (a) the full line 
is the prediction of the full 0{as) W+jet matrix element as implemented in Pythia. The 
dotted line is the prediction of the full DCM with /? = 2 and initial-state g qq splitting. 
In (h) the full line is as in (a) hut with the parton shower 0/ Pythia added after the first 
emission, the dashed line is Pythia using only parton showers and the dotted line is the 
same as in (a). 
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Figure 7: The average jet pseudo-rapidity rj vs. the rapidity of the W yw o,t the Tevatron. 
The jets are reconstructed with a cone algorithm using a radius of 0.7, and in each event a 
jet with E± > 20 GeV, and \r]\ < 3 on the oposite side in azimuth w.r.t. the W , is selected. 
In case of several such jets, the one with E± closest to the k± of the W is chosen. The full 
line is the leading order calculation as implemented in Pythia, the dashed line is the same, 
but with parton showers and fragmentation added, and the dotted line is the full DCM also 
with fragmentation added. (The "kinkiness" of the lines are due to limited statistics in the 
simulations.) 



Figure |^ is an attempt to reconstruct the measurement in ref. |T^0 on the generator level. 
As expected from eqs. (^) and (12), which are both symmetric around the W rapidity, 
the leading-order calculation gives a more or less linear correlation between the mean jet 
pseudo-rapidity and the rapidity of the W, although {rjj^t) 7^ Uq due to the smearing of the 
structure function convolution and the limited kinematical acceptance for jets. 

When parton showers and fragmentation are added to the leading-order calculation, the 
smearing is increased. Also, since the phase space available for emissions is larger on the 
side where the x of the incoming parton is smaller, the jets for large yw are "dragged" 
somewhat towards the centre, destroying the correlation. In the DCM, this dragging is 
more pronounced due to the ordering in the cascade as follows. 

In the parton shower in Pythia, each step in the backward evolution of the initial-state 
shower is ordered in both x and virtuality Q^; thus even if the phase space is larger on the 
small-x side of the W, the shower quickly runs out of phase space due to the ordering in 
(resulting also in an ordering in k"^) as in fig. The DCM, although ordered in k"^, is 
not ordered in x, or, if the final state partons are traced backwards in colour from the hard 
interaction, ordered in x but not in k"^ as in fig. ^o. In this respect, the DCM is similar to 
the BFKL evolution, and it gives a good description of the large transverse energy flows in 



^The details in tlie jet reconstruction may differ from that of ref. [ [ist . 
ambiguity in the yw determination is not taken into account here. 



In addition, the experimental 
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Figure 8: Example of paths, tracing emissions backwards from the hard interaction 
(yvyj Inm^j/) on the small-x side of high-rapidity W events for (a) the initial-state par- 
ton shower in Pythia, where the emissions are ordered both in x and k\, and (b) the 
DCM, where the emissions, although the cascade is ordered in k\, when traced backwards 
from the hard interaction in this way, are ordered in x but not in k'j_. 



small-x events at HERA, which has been suggested as a signal for the BFKL evolution |11]. 
Because of this, the DCM can better use the increased phase space on the small-x side of 
the W and the "dragging" effect is larger than for conventional parton showers in fig. 
and the result closer to, if not consistent with, the measurement in ref. |18]. 



The increase in transverse energy flow at small x found at HERA should also be visible 
at the Tevatron in high-rapidity W events. In fig. ^ the predictions for the E± flow from 
the parton shower model of Pythia^ and the DCM of Ariadne are shown for inclusive W 
events at the Tevatron for two W rapidity intervals. The two models are fairly similar at 
central W rapidities, while for high yw the CDM gives more transverse energy, despite the 
fact that the hard interaction scale {m^r ~ 6400 GeV^) is much larger here than at HERA 
((Q2) <ioo GeV2). 



In pp collisions we also have to worry about underlying events. In fig. ^, this does not give 
large effects since a large- i?_L jet is required, but for fig. ^ the underlying event would give 
an extra contribution to the E± flow. This extra contribution should however be evenly 
spread out in rj and independent of yw, and the differences between the parton shower 
and DCM approaches should survive. To take this contribution into account, the multiple 
interaction model implemented in Pythia |27| has been used. Note, however, that in the 
case of the DCM, only the qualitative features of the contribution are completely relevant, 
as the parameters of the multiple interaction model probably need to be retuned to fit the 
DCM. 



The differences between the parton shower and the DCM approaches are most significant 
when the E± flow is measured as a function of yw as in fig. 0, where the Ej_ flow, two 

®Since no high-fcx jets are required and the bulk of the events are at low k±w, the matrix element plus 
parton shower approach in fig. M is not adequate here. 
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(a) 0.0 < mi < 0.5 



(b) 2.0 < !j\y < 2.5 
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Figure 9: The transverse energy flow in inclusive W events at the Tevatron for (a) 0.0 < 
yw < 0.5 and (b) 2.0 < yw < 2.5. The full and dotted lines are the predictions of the DCM 
in Ariadne and of the parton shower in Pythia, respectively. 
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Figure 10: The transverse energy flow two units of rapidity "behind" the W for inclusive 
W events, at the Tevatron, i.e. for the interval 1.0 < yw < 1.5 the E± flow in the pseudo- 
rapidity interval —1.5 < rj < —1.0. The full line is the Pythia parton shower and the 
dotted line is the DCM without (a) and with (b) multiple interactions. 
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Figure 11: The transverse energy flow in the interval —2.5 < rj < —2.0 for W events at 
the Tevatron as a function of ywr- The full line is the Pythia parton shower and the 
dotted line is the DCM without (a) and with (h) multiple interactions. In (h) the dashed 
line is Pythia using 0{as) matrix elements with addition of parton showers and multiple 
interactions, and the dash- dotted line is the DCM with (3 = oo. 



units of rapidity away from the W , is shown for both models, with and without the multiple 



interaction model for the underlying event implemented in Pythia [27|. It is clear that 
the underlying event introduces an extra E_\_ for both models, but that the dependence of 
the E_\_ flow on yw is still different for the two models. As expected the Ex_ flow is more 
or less constant for the parton shower approach, but increases slightly for the DCM due to 
the increase in phase space at large yw- 

Similarly, looking at the E± flow in a fixed rapidity interval for varying yw, the DCM 
predicts a fairly constant value, while in the parton shower approach, the flow decreases 
with increasing yw as in fig. ITl|. To check that the differences are not due to the fact 



that the DCM has the correct 0{as) matrix element in the first emission, fig. lib also 
contains a line with the matrix element plus parton shower option in Pythia. It is clear 
that, although the E± flow is higher since only events with k±^w > 10 GeV are included, it 
has the same yw dependence as the plain parton shower approach. 



Also in fig. 11, the changing of the P parameter in the DCM is shown to have some effect 



on the E± flow; however, the dependence on yw is still the same. 



5 Conclusions 



The model presented in this paper is not perfect. The treatment of initial-state <?—>(?(? 
splitting is a bit foreign to the original Colour Dipole Model, and so is the transfer of 
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transverse recoil to the W for gluon emissions. But as a leading log approximation it should 
be just as good as the conventional parton shower approach, and it has the advantage of 
correctly describing also high-A;_|_ W production. In addition it gives the opportunity of 
studying effects of unordered parton evolution on the hadronic final state, and, although 
there are some uncertainties in the overall E± level due to multiple interactions and the (3 
parameter in the DCM, it would be very interesting to compare the predictions for the yw 
dependence of the E± flow presented in this paper with data from the Tevatron. 

Although only W production at the Tevatron has been discussed in this paper, the model 
can of course be applied to any Drell-Yan-like process in any hadron-hadron collision. And, 
as pointed out above, the initial-state g ^ qq splitting can be used also in deep inelastic 
lepton-hadron scattering. 

The model presented here is also the "last piece" to complete the DCM description of QCD 

cascades for all standard processes in e+e", ep and pp collisions. This is reflected in the 
fact that Ariadne now is fully interfaced to all hard sub-processes in Pythia. However, 
some care must be taken when using the multiple interaction model of Pythia together 
with the DCM, as mentioned above. 
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